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Abstract 

Freeze out of particles across three dimensional space-time hypersurface is discussed 
in a simple kinetic model. The final momentum distribution of emitted particles, 
for freeze out surfaces with space-like normal, shows a non-exponential transverse 
momentum spectrum. The slope parameter of the pt distribution increases with 
increasing pt, in agreement with recently measured SPS pion and h~ spectra. 
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1 Introduction 

In continuum and fluid dynamical models, particles, which leave the system 
and reach the detectors, can be taken into account via freeze out (FO) or final 
break-up schemes, where the frozen out particles are formed on a 3-dimensional 
hypersurface in space-time. Such FO descriptions are important ingredients of 
evaluations of two-particle correlation data, transverse-, longitudinal-, radial- 
, and cylindrical- flow analyses, transverse momentum and transverse mass 
spectra and many other observables. The FO on a hypersurface is a disconti- 
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nuity where the pre-FO equilibrated and interacting matter abruptly changes 
to non-interacting particles, showing an ideal gas type of behavior. 



The general theory of discontinuities in relativistic flow was not worked out 
for a long time, and the 1948 work of A. Taub (|J) discussed discontinuities 
across propagating hypersurfaces only (which have a space-like unit normal 
vector, dd^db^ = — 1). Events happening on a propagating, (2 dimensional) 
surface belong to this category. An overall sudden change in a finite volume 
is represented by a hypersurface with a time- like normal, da^da^ = +1. The 
freeze out surface is frequently a surface with time like normal. In 1987 Taub's 
approach was generalized to both types of surfaces (0), making it possible to 
take into account conservation laws exactly across any surface of discontinuity 
in relativistic flow. When the EoS is different on the two sides of the freeze 
out front these conservation laws yield changing temperature, density, flow 
velocity across the front (|3|; f|; |5|; |6|; |7|). 



2 Conservation laws across idealized freeze out fronts 

The freeze out surface is an idealization of a layer of finite thickness (of the 
order of a mean free path or collision time) where the frozen-out particles are 
formed and the interactions in the matter become negligible. 

To use well-known Cooper-Frye formula (|||) 



we have to know the post-FO distribution of frozen out particles, fFo( x , P] T, n, u u ), 
which is not known from the fluid dynamical model. To evaluate measurables 
we have to know the correct parameters of the matter after the FO discon- 
tinuity! The post freeze out distribution need not be a thermal distribution! 
In fact jpo should contain only particles which cross the FO-front outwards, 
p^da^ > 0, so if da^ is space- like this seriously constrains the shape of jpo- 
This problem was recognized in recent years, and the first suggestions for the 
solution were published recently (U |5| || 

If we know the pre freeze out baryon current and energy-momentum tensor, 
N{f and Tq U , we can calculate locally, across a surface element of normal 
vector d&L 1 the post freeze out quantities, iV M and T^, from the relations 
(Q; |): [iV M da^} = and [T^ da^} = 0, where [A] = A — A . In numerical 
calculations the local freeze out surface can be determined most accurately 
via self-consistent iteration (f|; |9]). 
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3 Freeze out distribution from kinetic theory 

We present a kinetic model simplified to the limit where we can obtain a post 
FO particle momentum distribution. Let us assume an infinitely long tube 
with its left half (x < 0) filled with nuclear mater and in the right vacuum is 
maintained. We can remove the dividing wall at t — 0, and then the matter 
will expand into the vacuum. By continuously removing particles at the right 
end of the tube and supplying particles on the left end, we can establish a 
stationary flow in the tube, where the particles will gradually freeze out in an 
exponential rarefraction wave propagating to the left in the matter. We can 
move with this front, so that we describe it from the reference frame of the 
front (RFF). 

We can describe the freeze out kinetics on the r.h.s. of the tube assuming 
that we have two components of our momentum distribution, ff ree (x,p) and 
fint(x,p). However, we assume that at x = 0, ff ree vanishes exactly and fi nt 
is an ideal Jiittner distribution, then f int gradually disappears and ff ree grad- 
ually builds up. 

Rescattering within the interacting component will lead to re-thermalization 
and re-equilibration of this component. Thus, the evolution of the component, 
fint is determined by drain terms and the re-equilibration. 

We use the relaxation time approximation to simplify the description of the 
dynamics. Then the two components of the momentum distribution develop 
according to the coupled differential equations: 

d x f int (x,p)dx = -&(p' l da fl )^ £ fi n t(x,p)dx+ 

(2) 

+ [feq(x,p) - f int (x,p)] jjdx, 

d x f free (x,p)dx = +®(p^dcr^)^£f int (x,p)dx. (3) 

Here cos^ = p x /p in the RFF frame. The first (loss) term in eq. (2) is an 
overly simplified approximation to the model presented in ref. (||). It expresses 
the fact that particles with momenta orthogonal to the FO surface (cos 6^ = 1) 
leave the system with bigger probability than particles emitted at an angle. 
The interacting component of the momentum distribution, described by eq. 
(2), shows the tendency to approach an equilibrated distribution with a relax- 
ation length A'. Of course, due to the energy, momentum and particle drain, 
this distribution, f eq (x,p) is not the same as the initial Jiittner distribution, 
but its parameters, n eq (x), T eq (x) and [x\ change as required by the con- 
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servation laws. 



In this case the change of the conserved quantities caused by the particle 
transfer from component int to component free can be obtained in terms of 
the distribution functions as: 

dN » = _dx I <$v ^ cosB ( 
A J Po 

and 

= _dx I fP p , fQ ^ d&) CQs6 f ( (5) 
A J po 

Due to the collision or relaxation terms T^ u and change, and this should 
be considered in the modified distribution function f int (x,p). 



3.1 Immediate re-thermalization limit 

Let us assume that A' <C A, i.e. re-thermalization is much faster than particles 
freezing out, or much faster than parameters, n eq (x), T eq (x) and u% q (x) change. 
Then f in t(x,p) » f eq (x,p), for A'«A. 

For f eq {x, p) we assume the spherical Jiittner form at any x including both pos- 
itive and negative momentum parts with parameters n(x), T(x) and u RFG (x). 
(Here u RFG (x) is the actual flow velocity of the interacting, Jiittner compo- 
nent, i.e. the velocity of the Rest Frame of the Gas (RFG) (|])). 

In this case the change of conserved quantities due to particle drain or transfer 
can be evaluated for an infinitesimal dx. The changes of the conserved particle 
currents and energy-momentum tensor in the RFF, eqs. (4, 5) are given in ref. 
(H). The new parameters of distribution f int , after moving to the right by dx 
can be obtained from dN- 1 and dTf v . The differential equation describing the 
change of the proper particle density is (0): 

dui(x) = uf jRFG (x) dNi^(x) . (6) 

Although this covariant equation is valid in any frame, dN- 1 are calculated in 
the RFF §). 

For the re-thermalized interacting component the change of Eckart's flow ve- 
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locity is given by 



ni{x) 



where Af^(x) = — u^ RFG (x) RFG (x) is a projector to the plane orthog- 
onal to uf RFG (x), while the change of Landau's flow velocity is (^) 



e* + Pi 



au i,L,RFG\ x ) — _ , n • v°/ 



Although, for the spherical Juttner distribution the Landau and Eckart flow 
velocities are the same, the change of this flow velocity calculated from the loss 
of baryon current and from the loss of energy current are different duf ERFG (x) ^ 
du^ LRFG (x) . This is a clear consequence of the asymmetry caused by the 
freeze out process as it was discussed in ref. (H), i.e., the cut by Qip^da^) 
changes the particle flow and energy-momentum flow differently. This prob- 
lem does not occur for the freeze out of baryonfree plasma, and we have only 
du? L . 
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Fig. 1. The local transverse momentum (herepa;) distribution for baryon free, mass- 
less gas at p y = 0, x = 100A and To = 130 MeV. The transverse momentum spec- 
trum is obviously curved due to the freeze out process, particularly for large initial 
flow velocities. The apparent slope parameter increases with increasing transverse 
momentum. This behavior agrees with observed pion transverse mass spectra at 



SPS fllj; y) 



The last task is to determine the change of the temperature parameter of fi nt . 
From the relation e = u^T^Up we readily obtain the expression for the change 
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of energy density 

dei{x) = u^RFaix) dT^ u (x) u U)i)RFG (x) , (9) 



and from the relation between the energy density and the temperature (see 
Chapter 3 in ref. (§)), we can obtain the new temperature at x + dx. Fixing 
these parameters we fully determined the spherical Juttner approximation for 

fint- 

The application of this model to the baryonfree and massless gas gives the 
following coupled set of equations: 



d\nT 

dx 



du^ 
dx 



(10) 



-u*u T av ] u v 



Here we use the EoS, e = er^T 4 , the definition dT^ v = —dxr^T 4 , and x 
is measured in units of A. 

Now we can find the distribution function for the noninteracting, frozen out 
part of particles according to equation (3). The results are shown in Fig. 1. 
We would like to note that now f int {x,f>) does not tend to the cut Juttner 
distribution in the limit x — > oo. Furthermore, we obtain that T — > 0, when 
x — > oo (I). So, fmt{x,p) = (2^)3 exp[(/x - p v u v ) /T] -> 0, when x -> oo. Thus, 
all particles freeze out in the present model, but such a physical FO requires 
infinite distance (or time). This second problem may also be removed by using 
volume emission model discussed in (|7|). 



4 Conclusions 



In a simple kinetic model we evaluated the freeze out distribution, ff ree (x,p), 
for 

stationary freeze out across a surface with space-like normal vector, da^da^ < 
0. In this model particles penetrating the surface outwards were allowed to 
freeze out with a probability ~ cos 9p, and the remaining interacting compo- 
nent is assumed to be instantly re-thermalized. The three parameters of the 
interacting component, f int , are obtained in each time step. The density of 
the interacting component gradually decreases and disappears, the flow ve- 
locity also decreases and the energy density decreases. The temperature, as 
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a consequence of the gradual change in the emission mechanism, gradually 
decreases. 

The arising post freeze out distribution, fj ree is a superposition of cut Jiittner 
type of components, from a series of gradually slowing down Jiittner distri- 
butions. This leads to a final momentum distribution, with a more dominant 
peak at zero momentum and a forward halo, Fig. 1. In this rough model a 
large fraction (~ 95%) of the matter is frozen out by x = 3A, thus, the dis- 
tribution ff ree at this distance can be considered as a first estimation of the 
post freeze out distribution. One should also keep in mind that the model 
presented here does not have realistic behavior in the limit x — > oo, due to 
its one dimensional character. 



These studies indicate that more attention should be paid to the final freeze 
out process, because a realistic freeze out description may lead to large p t en- 
hancement ( |10t |TTD as the considerations above indicate (Fig. 1). For accurate 
estimates more realistic models should be used. In case of rapid hadronization 
of QGP and simultaneous freeze out, the idealization of a freeze out hypersur- 
face may be justified, however, an accurate determination of the post freeze 
out hadron momentum distribution would require a nontrivial dynamical cal- 
culation. 
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